Chiral condensate in hadronic matter 
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The finite temperature chiral condensate for 2 + 1 quark flavors is considered in the framework of 
the hadron resonance gas model. This requires some dynamical information, for which two models 
are employed: one based on the quark structure of hadrons combined with the Nambu-Jona-Lasinio 
approach to chiral symmetry breaking, and one originating from gauge/gravity duality. Using these 
insights, hadronic sigma terms are discussed in the context of recent first principles results following 
from lattice QCD and chiral perturbation theory. For the condensate, in generic agreement with 
lattice data it is found that chiral symmetry restoration in the strange quark sector takes place at 
higher temperatures than in the light quark sector. The importance of this result for a recently 
proposed dynamical model of hadronic freeze-out is outlined. 
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I. INTRODUCTION 

Spontaneous chiral symmetry breaking is, apart from 
color confinement, the most important physical aspect 
of strong interactions. The fact that one observes mass 
splittings of chiral partners in the hadron spectrum 
and that pions have properties attributed to Goldstone 
bosons strongly suggests that chiral symmetry is sponta- 
neously broken in the vacuum. These, and other theoret- 
ical arguments [l[ , imply that in the vacuum there exists 
chiral condensate 1 giving rise to an expectation value of 
the bilinear fermionic operator ipip. Dynamical details of 
this phenomenon, which is inherently non-perturbative in 
nature, are part of the long standing problem of strong 
interactions, but in the course of time different model 
mechanisms for all its different aspects have been devel- 
oped. 

As temperature and/or baryon density is increased, 
thermal hadron excitations, because of their quark sub- 
structure, will affect the vacuum condensate causing its 
melting and eventually vanishing at the transition line to 
the chirally symmetric phase. Microscopic quantification 
of this phenomenon comes from first principles lattice 
QCD (1QCD) simulations and confirms the intuitive pre- 
dictions Q. 

To get physical insight into this effect for low temper- 
atures (and densities) one can use the hadron resonance 
gas (HRG) model @, which was previously successfully 
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1 Note however, that it may be possible to reinterpret the evidence 
in the spirit of 0,0] where all effects of chiral symmetry breaking 
are attributed to hadronic states rather than to the QCD ground 
state. 



applied to give physical interpretation of QCD data (HQ 
as well as a description of the abundances of particles 
produced in heavy ion collisions at very different center 
of mass energies @, Q in terms of freeze-out parameters. 
The assumption underlying this approach is that for con- 
ditions below the QCD transition line the system is com- 
posed of non-interacting hadronic degrees of freedom and 
so the partition function is that of an ideal mixture of free 
quantum gases. To have a reliable physical description 
of the system one needs to take into account all hadron 
resonances with masses up to ~ 2 GeV. 

To calculate the condensate in this framework it is 
necessary to know the dependence of hadron masses on 
the current quark masses. This, apart from the Gold- 
stone boson octet, is not straightforward to determine 
and either requires some assumptions about the under- 
lying dynamics or is the result of a phenomenological 
fit. Approaches, which can be regarded as based on first 
principles, are chiral perturbation theory (ChPT) lat- 
tice QCD simulations [l(| and Dyson-Schwinger equa- 
tions (DSE) [0, HI]- They provide a consistent picture 
of hadrons with a reliable account of the quark mass de- 
pendence. However, in the ChPT framework there are 
still large uncertainties concerning for example the nu- 
cleon strange sigma term fl3| for which, when different 
orders of app roximation are considered, even the sign is 
not clear |14| . 

This article explores the consequences of various 
hadron mass formulae proposed recently and compares 
them with the results mentioned above. 

One set of mass formulae which were used in the analy- 
sis reported here comes from a new model based on quark 
counting and is a generalization of what was proposed by 
Leupold [IH a few years ago. In Leupold's scheme hadron 
masses were assumed to be linear in the current quark 
masses. This approach was used (and generalized) in fl6| . 
In the present work a further step is taken: it is assumed 
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that the dependence of hadron masses on the current 
quark mass arises solely due to the dependence of con- 
stituent masses of valence quarks. The response of these 
constituent masses to the change in the current quark 
mass is determined based on the Nambu-Jona-Lasinio 
(NJL) model ftft. In this way a fairly good description 
of the hadronic sigma terms is obtained. The only flaw 
is that the sea-quark contributions arc neglected entirely, 
which, for example, leads to the vanishing of the nucleon 
strange sigma term. 

The second approach considered in this paper is the 
use of baryon mass formulae which were obtained in a 
large N |18| holographic model of QCD due to Sakai and 
Sugimoto |19| . The last ten years have witnessed a lot of 
progress coming from gauge/gravity duality allowing for 
valuable insights into the dynamics of strongly coupled 
gauge theories. Recent developments have made it possi- 
ble to study in a quasi-analytical way theories which have 
very realistic properties. The spectrum of mesons and 
chiral symmetry breaking in the chiral limit was stud- 
ied in [19j and baryon masses, static properties (such as 
magnetic moments or charge radii) [20l . |21| were found 
to be in qualitative agreement with experiment. Also 
form factors 2l| agree quite well with the data. Further 
progress was made with the extension to finite current 
quark masses [22| (see [23| for an alternative construc- 
tion) where for example Gell-Mann-Oakes-Renner rela- 
tions for the pseudoscalar octet have been demonstrated. 
The impact of finite current quark masses on the spec- 
trum of nucleon octet and delta decuplet baryons has 
been considered in the two flavor case [24| and for 2 + 1 
flavors [25j j with nontrivial results. The leading order cor- 
rections are proportional to the squares of Goldstone bo- 
son masses and determined in a similar way as the leading 
order of ChPT (l4| . The results are in good agreement 
with the data and other theoretical expectations in the 
light quark sector, while the contribution of the strange 
quark is overestimated by the model. It is very likely 
that going beyond the leading order in the expansion in 
powers of the current quark masses will give more rea- 
sonable results (as is the case in ChPT). Also, at leading 
order, vector mesons were argued not to receive mass cor- 
rections from finite current quark mass [24j . The mass 
formulae obtained in [24|, HH make it possible to esti- 
mate sigma terms, including those for the nucleon octet 
and delta decuplet. This is then used to calculate the 
chiral condensate in the framework of the HRG model 
and the results are compared with calculations based on 
chiral perturbation theory. 

In the context of DSE studies [IJ IH1 sigma terms for 
the two light quark flavours have been considered. In 
addition to the nucleon and delta baryons also vector 
mesons were included. Due to the p — tttt and p — ir — lu 
one gets sigma term of the p-meson. The w-meson has 
no pion loop dressing and therefore only lu — pir coupling 
remains. Because in the DSE approach strange sigma 
terms were not included we do not use this results as 
a base for the calculation of the chiral condensates of 



interest. We will only use it as a reference point to other 
calculations. 

The importance of the hadronic contribution to the 
melting of the chiral condensate was appreciated in a 
model for the freeze-out stage of heavy-ion collisions, 
where it was related to the Mott- Anderson dclocalisa- 
tion of hadrons [l(|. The model is based on several 
detailed assumptions on the evolution of the matter af- 
ter heavy ion collision and hadron interactions, but the 
main point is that freeze-out phenomena are assumed 
to take place in the hadronic phase and are entirely at- 
tributed to the hadron dynamics. In general, each hadron 
is assigned a medium dependent radius r^(T, /is), which 
is then related in an universal way to the chiral con- 
densate. Hadron-hadron reactions are described by the 
Povh-Hufner law [26j and in consequence the cross sec- 
tion is determined by the medium dependence of the 
condensate. As the temperature decreases, the mean 
time between the interactions is getting larger, since it is 
inversely proportional to the reaction cross-section and 
hadron density (in the relaxation time approximation). 
At some point this becomes larger than the rate of ex- 
pansion and reactions between hadrons stop to change 
the final composition. The freeze-out parameters p B and 
are determined by the equality of both time scales. 
However in Rcf. fl6j only the light quark condensate was 
considered, so one of the possible improvements of the 
model is to include also the strange sector. This is one 
of the motivations for the present studies. 

The organization of the paper is as follows: in section 
HT1 the generic theoretical setup is described and the rel- 
evant quantities used in further calculations are defined. 
This section also reviews some thermodynamic quanti- 
ties computed in the HRG model and highlights very 
good agreement with lattice computations. These consid- 
erations do not require any detailed assumptions about 
hadron dynamics. On the other hand, the computation 
of the chiral condensate strongly depends on hadron mass 
formulae expressed in terms of current quark masses as 
discussed in section [TTT1 This dependence is captured by 
the hadronic sigma terms. In section IPV1 we describe our 
baseline which are results obtained within ChPT as the 
low energy effective theory of QCD. In section [V] mass 
formulae based on constituent quark structure are pre- 
sented and contrasted with previously established para- 
metric dependence [13, [28[ and with first principles re- 
sults. Section IVII contains novel results following from 
the Sakai-Sugimoto holographic model together with a 
discussion in the light of lowest order ChPT results. Sec- 
tion [VTT] contains conclusions, some discussion and pos- 
sible open directions. 



II. HADRON RESONANCE GAS MODEL 

_ The hadron resonance gas model implements the idea 
[j| that QCD thermodynamics in the hadronic phase can 
be described as a multicomponcnt ideal hadron gas. For 
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very low temperatures the dominant degrees of freedom 
are pions and kaons and due to the Goldstone theorem 
their interactions are weak. Therefore in the first approx- 
imation they can be considered as free particles. As the 
temperature and/or density is increased, contributions 
from heavier hadrons become important. Because strong 
interactions are of finite range in the thermodynamical 
limit of infinite volume V — > oo the grand canonical po- 
tential 2 can be expressed as a sum of contributions from 
free hadrons [301 



n(T, {&}) = n + sihrg(t, {fti}) 



(1) 



In the above formula £lo is the vacuum part, whose de- 
tailed form is irrelevant for the following considerations, 
and the medium dependent part contains contributions 
from mesons and baryons 

n HRG (T, fa}) = Q M (T, fa}) + Q B (T, fa}) . (2) 

Here {fJ-i} is the set of chemical potentials corresponding 
to conserved charges such as baryon number B, electric 
charge Q, isospin / 3 and strangeness S. The free meson 
contribution reads 



-z M e 



~ 0Em ), (3) 



r Ail. 



while the free baryon contribution is 

O fl (T,{ Mi }) = -^d fl J j^Tln(l + z B e-^) , 

B _ (4) 

where Ei = y k 2 + mf and ds and d,M count the degen- 
eracy of hadrons. Fugacitics arc defined by 



(5) 



where the index a runs over all conserved charges in the 
system, X a is the corresponding charge and (3 = l/T is 
the inverse temperature. Although inclusion of chemical 
potentials is straightforward in the HRG approach, for 
much of this paper they are all set to zero. The residual 
repulsive interactions can be taken into account , e.g ., by 
the van der Waals excluded volume corrections 
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All thermodynamic quantities, such as equations of 
state for pressure and energy density as well as mate- 
rial properties such as the speed of sound, can be ob- 
tained from the grand canonical thermodynamical po- 
tential f2(T, {fJ-i}). The pressure is given by 



p = -fi(T, fa = 0}) 



(6) 



2 Since only homogeneous systems are considered here, the symbol 
Q denotes the grand canonical potential as usually defined in 
statistical mechanics divided by the volume. 
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FIG. 1: Energy density and pressure for HRG and lattice data 



and the energy density is 



£ = 



d\pn(T, {^ = o})] 

0(3 



(7) 



Fig. Q] shows the equations of state as obtained for the 
HRG and compares it with recent 1QCD simulations [HJ 
normalized to the T 4 (Stefan-Boltzmann) behaviour of a 
massless ideal gas. There is clearly an excellent agree- 
ment for temperatures up to ~ 170 MeV which means 
that the dominant effect in that range of temperatures 
comes from the excitation of hadronic degrees of freedom 
rather than from their interactions. This is a very well 
known effect d, Q . Agreement for higher temperatures 
can be obtained when medium modifications of hadronic 
states are taken into account. For example in [33| it was 
demonstrated that inclusion of state dependent hadronic 
width Th(T) taken on the inverse collision time of the 
Mott- Anderson freeze-out (lrj | and proper introduction of 
quark-gluon degrees of freedom based on the Polyakov- 
loop extended Nambu-Jona-Lasinio (PNJL) model nicely 
reproduces lattice QCD data in the whole temperature 
range. 

The velocity of sound is given by 



c 2 = 



dp e + p ( de 

It 



de 



T 



(8) 



where the second equality holds only for zero chemical 
potentials. Its temperature dependence is shown in Fig. [2] 
for the HRG model compared to 1QCD data [32j |. 

Qualitatively in the HRG model the increase for low 
temperatures is related to the appearance of a large num- 
ber of light degrees of freedom. When heavier hadrons 
are exited, they contribute considerably to the energy 
density but almost nothing to the pressure, which leads 
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FIG. 2: Sound velocity for HRG and lattice data [13] 



The derivatives taken with respect to the current quark 
masses and lead to the generic formulae 



<qq> = <qq> +y2^n M (T,{^}) (12) 



<ss> = < ss > + Y\— n M {T,{m}) (13) 



where the scalar densities of mesons and baryons have 
been introduced as 



n M (T,{fii}) = -5§ / dkk 2 ^- l, - 



(14) 



to the characteristic dip. For high temperatures, be- 
cause the number of states included is finite, there is an 
approximately constant behavior approaching the mass- 
less gas limit c 2 . = 1/3 only for very high temperatures. 
On the other hand, for 1QCD the dip is an indicator of 
the crossover transition where the sound velocity should 
be strictly zero. For high temperatures lattice data ap- 
proach the masslcss limit, which is consistent with the in- 
terpretation of deconfinement to a masslcss quark-gluon 
medium. 

The importance of the speed of sound for the phe- 
nomenology of heavy ion collisions was noticed by 
Florkowski et at. [H in the context of the HBT 
puzzle. 



III. CHIRAL CONDENSATE AND SIGMA 
TERMS 



Using the standard formula 

9fi(T,{M) 



< qq >- 



dmn 



(9) 



one obtains the quark-antiquark condensate in the light 
and strange flavor sector respectively 



< qq >=< qq > + 



(10) 



n B (T, {m}) = 



2tt 2 



dkk 



2 m B 



E B ZgV^ 



1 



(15) 



and the response of hadron masses to changes in the cur- 
rent quark mass of flavor / = u 1 d 1 s, gjy. is captured 
by the hadron sigma terms 



dm h 
dmj 



(16) 



Thus, for every hadron state, there arc different sigma 
terms related to contributions from quark flavors consti- 
tuting the hadron. 

The above formulas are valid for the non-interacting 
gas. In the light flavour sector effects of meson-meson 
and pion-nucleon interactions as described by ChPT were 
implemented in |36j ] . 

In the following, isospin symmetry is assumed, setting 
the light quark masses m q = m u = rj 5.5 MeV 
and the light quark condensate < uu >o=< dd >o=< 
qq > = (-240) 3 MeV 3 . Analysis based on 1QCD and 
QCD sum rules together with the low energy theorem 
for the correlation functions allows one to estimate the 
ratio of strange to light quark condensates to be 0.8 ±0.3 
[37l | (other estimates give 0.75 ± 0.12 [38| but note also 
recent explicit lattice calculations (HI). One can under- 
stand this hierarchy of condensates using the spectral 
representation of the expectation value for the quark of 
current mass m/ [H, [40l ] 



-2m 



dX- 



P(A) 



A 2 



(17) 



< ss >=< ss > H = • (11) 

om s 



and noting that the spectral integral is increasingly sup- 
pressed with the higher current quark mass, thus lower- 
ing the value of the quark condensate. Furthermore, the 
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characteristic length scale related to the quark-antiquark 
condensate can be taken as l/rrif which is smaller for 
greater masses. This implies that the medium effect - 
expressed as screening length - will affect heavier quark 
condensates at higher temperatures. This can also be un- 
derstood as arising from the fact that the contribution to 
the strange quark condensate - and its melting - comes 
from strange hadrons, which are fewer in number than 
hadrons containing light quarks. 

Another important quantity, an approximate order 
parameter for the deconfinement phase transition, is 
Polyakov loop. It is very well studied in 1QCD and re- 
cently it has been addressed within the HRG framework 
[1H . Good agreement with the lattice data was found in 
the temperature range 150 MeV< T < 190 MeV. 



Taking the derivative of the above equations with re- 
spect to the light quark masses one finds 



dm 2 
dm a 



< qq >o 



(20) 

and similarly for the derivatives of the kaon mass with 
respect to m q (and m s ): 



dm a , s 



< qq > + < ss > 



2/£(l-2^ 

< qq >p + < ss >p 
2/^ 



(21) 



1 + 2k- 



K 



f-2 

J 7T 



IV. HADRON MASSES IN CHIRAL 
PERTURBATION THEORY 

As explained above, the finite temperature behaviour 
of chiral condensates in the HRG approach is determined 
by the sigma terms, which express the dependence of 
hadron masses on the current quark masses. This paper 
explores two models for this dependence (discussed in the 
following sections), and compares them with the results 
of chiral perturbation theory, which is the subject of this 
section. 

It is natural to start with the Goldstone bosons octet. 
In this case the relevant mass formula is the Gell-Mann- 
Oakes-Renner (GMOR) relation, which takes the form 

m 



r2 2 



f 2 

J 7T 



1 - K-jf = - < qq >o (m u + m d ) , (18) 



!k™k 1 



< K 



f 2 



< qq > + < ss > 



(m q +m s ) . 
(19) 

These formulae include next to leading order corrections 
expressed in terms of the parameter k = 0.021 ± 0.008 
[37j . If one assumes < ss > = 0.8 < qq >o, = 
92.4 MeV, f K = 113 MeV (which gives f K /h ~ 
1.22 0) and m s = 138 MeV, then one finds (m q + 
m s )/m s fa 1.040 as compared to the lattice choice Q 
(m, q + m s )/m s w 1.036. 



In ChPT mass formulae for the ground state baryons 
can be also computed in the Nf = 2 j44j and Nf =2 + 1 



cases 



141. At lowest order the shift due to the finite 



current quark mass is proportional to the square of the 
Goldstone boson mass [lj]. The lowest order contribu- 
tions to the baryon masses read 



Mb = Mi - 



f>=TT,K 



(22) 



where £b,0 are expressed in terms of the parameters of 
the low energy Lagrangian. Sigma terms following from 
this formula are listed in the table |U In section PVTl these 
results will compared with the holographic mass formu- 
lae. 

Higher order contributions are given by the loop cor- 
rections to baryon self energies and are evaluated using 
different rcgularization schemes. All of them are carefully 
compared to the recent lattice data. It is interesting to 
note that the strange nucleon sigma term turns out to be 
negative at next-to-leading order (NLO) cr^ « —4 MeV 
[l4| which means that the nucleon mass decreases if the 
strange quark mass is raised. On the other hand differ- 
ent ChPT studies give at N 3 LO erf w 130 MeVfH which 
means that higher order corrections can be relevant and 
the existing answers must be regarded as somewhat ten- 
tative. Firs principle lattice simulations with Nf =2 + 1 
dynamical quark flavours give trf = 49 + 25 MeV [45j] . 





N 


A 


E 




A 


E* 




Q~ 


(T.s 


36.2527 
162.474 


19.8647 
437.693 


10.7535 
590.705 


27.8369 
317.096 


-1.07895 
789.418 


15.5728 
523.058 


3.30878 
729.019 


-8.95528 
934.981 



TABLE I: Sigma terms for the lowest lying baryons in the leading order ChPT in MeV 
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In this context an interesting quantity to look at is 
strangness content of the nucleon which was estimated 
in the form Q 



the hidden strange mesons. The strangess counting pa- 



2 < p|ss|p > 
< p\uu + dd\p > 
,2 / i 

n - 1 



(23) 



&ttN 



>K 



2 W ' 



dm 



N 



0.21 



where \p > is a nucleon state of momentum p. This is 
similar to the famous Wroblewski factor [46| introduced 
in HIC and pp collisions for the quantification of strang- 
ness production. For the SU(3) symmetric case y = 1 
while when there are no strange quark pairs y = 0. Sec- 
ond equality of the Eq. [23] is quite generic and relies only 
on the Hellman-Fyenman theorem and tree level GMOR 
relations. Its importance lies in the fact that making 
some statement about the nucleon strange sigma term 
is in fact equivalent to making a statement about the 
strangness contribution to the nucleon. 



V. CONSTITUENT QUARK PICTURE 

The model described in this section is based on the va- 
lence quark structure of hadrons and is a nontrivial gen- 
eralization of the formulae used in [l6[ for the light quark 
condensate (following earlier work by Leupoldjla]). This 
model is also compared with a another approach which 
gives a parametric dependence of hadron masses on the 
pion mass [13, [28| and was previously used for the calcu- 
lation of both light and strange quark condensates [2!| . 
In these two models mass formulae arc in principle given 
for all the hadron states and so the sums over mesons 
(M) and baryons (B) which appear in the HRG model 
take into account all states up to mass ~ 2 GeV. 

The scenario introduced here assumes that baryon and 
meson masses scale as 

m B = (3 - N s )M q + N S M S + k b , (24) 



m M = (2 - N s )M q + N S M S + K M ■ (25) 

The quark masses in these mass formulae are the dynami- 
cal (constituent) ones and are denoted by M q for the light 
quarks and by M s for the strange quark. The parameter 
N a measures the strangess content of the hadron and the 
quantities Kb, k-m depend on the state, but not on the 
current quark masses. For the open strange hadrons N s 
is simply the number of strange (anti-strange) quarks. 
For hidden strange mesons - such as for example the rj 
or hi state - it is modified by the square modulus of the 
coefficient of the ss contribution to the meson wave func- 
tion. There are two possible wave function assignments 
related to the flavor singlet ipo = -^={uu + dd + ss) and 

flavor octet tpg = -^=(uu + dd — 2ss) wave functions for 



r (0) 



2/3 for singlet and vi 8) =4/3 for the 



ramcter iV s 
octet. 

An easy retrospection shows that the baryon octet 
Gell-Mann-Okubo relations: 3Af A + M s = 2(M N + M s ) 
is translated into constraint on the state dependent con- 
tributions: 3ka + ks = 2(kn + ks) 

Two further simplifying assumptions arc made: ex- 
cited states are assumed to have the same wave functions 
as their respective ground states, and any possible mixing 
between octet and singlet states (such as r\ — rf mixing) 
is neglected. 

The dynamical (constituent) quark masses M q and M s 
appearing in Eq. (|24 p .(|25 p are a way of partially ac- 
counting for the dynamics of strong interactions. For the 
purposes of computing the condensates only the depen- 
dence of these constituent masses on the current quark 
masses is relevant. This dependence is taken from NJL 
model, where dynamically generated mass changes by 
AM q = 12.5 MeV as the quark mass is turned on from 
zero in the chiral limit to m q = 5.5 MeV. This gives the 
nucleon sigma term cr/v = 37.5 MeV. For the strange 
quark mass the value of the dynamical quark mass is 
M s = 587.4 MeV for m s = 140.7 MeV which gives 
AM, = 227.4 MeV |47|. This valence quark counting 



implies that the strange contribution to the nucleon is 
zero which is an approximation hard to control. For the 
A baryon which has one strange quark the same argu- 
ments as above lead to the estimate = 252.9 MeV. 
The resulting sigma terms are shown in Figs. 13141 

In principle the scaling of Eqs. (|24p. (|25l) will be cor- 
rected by various effects such as contributions of the sea 
quarks, which one would expect to give a logarithmic 
correction ln(m g / J^-qcd) 2 at one loop order. 

At this point we would like to discuss a previously used 
parametrization for the dependence of hadron mass on 
the explicit breaking of chiral symmetry, i.e., on the pion 
mass squared |27H29l | 

drrih Ah , , 

(26) 



dmi. 



m h 



where Ah ~ 0.9 — 1.2. This parametrization is particu- 
larly interesting when using lattice QCD data for hadron 
masses which are determined at different quark masses 
corresponding to squared pseudoscalar meson masses 
(to^.). The state-dependent coefficient Ah can be used 
to replace the sigma terms in the condensate formulas 
([IB"]). dHJ) according to (assuming k = 0) 



= ^A h 
m h 



(27) 



for the light quark sigma-terms. Below we will use also 
this formula to translate sigma terms calculated within 
CQP to estimate A h . Using the GMOR relation (fTS]). 
one can write for the light quark condensate in a HRG 
the compact expression 



<qq> = < qq >o 1 



A av ntot 



(28) 
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where the averaged Ah coefficient is introduced as 
. J2h={Ai},{B} Ahrih/mh 



J2h={M},{B} n h/m h 



while 



TO re d 



h={M},{B} 



n h /m h 



Ylh={M},{B} n h 



(29) 



(30) 



is the weighted reduced mass and n to t = J2h={M} {b} n h 
the total scalar density of hadrons. Note, that A av and 
?7i rc d are temperature dependent. 

The Eq. (f2"5| gives a compact expression for the modi- 
fication of the light quark condensate in a HRG medium. 
Since n to t and m r od are model independent character- 
istics of the HRG, the evaluation of the medium depen- 
dence requires solely the determination of A av for a given 
model. 

The reduced mass defined in equation ([30)) is analogous 
to the reduced mass fi rc d used in many particle system. 
The latter obeys two inequalities mughtest/'i < A*red < 
^lightest, where miightcst is the lightest mass in the system 
of n particles. Those inequalities have direct analogy in 
our case and read < m re d < m 7r n to t(r)/n 7r (T) where 
pion is the lightest hadron and n v {T) is scalar density of 
pion. 

We exemplify this for the simple quark counting model 
with the mass formulas of Eqs. (|2"4")l and (|2~5j) for which 
we have already given the sigma terms. The correspond- 
ing values of the Ah coefficient as a function of hadron 
mass are shown in Fig. [5] For this model, the averaged 
value (|2"9")l comes out to be temperature dependent and 
its behaviour is shown in Fig. [S] for three different upper 
limits in the included hadron spectrum. 

The straight line structures of Fig. [5] reflect the fact 
that different hadrons admit different flavour structure 
and the assumption that excited states have the same 
structure as their respective ground states. 

Figure [7] shows chiral condensates calculated with the 
two mass formulae described above. What is apparent 
is that in the quark counting scenario there is a more 
pronounced difference between the light and the strange 
condensates. In I29f it was found that for the paramet- 
ric mass formulae [27t [28j at the temperature where the 
light condensate vanishes the strange condensate is 0.4 
of its vacuum value. The temperature where the light 
condensate vanishes is about T ~ 190 McV. In con- 
trast for quark counting scheme used here this ratio is 
~ 0.78. The temperature where the light quark 
vanishes is T « 168 McV. This difference comes from 
the fact that taking into account sea quark effects dimin- 
ishes the difference between contributions from strange 
and non-strange hadrons. For example nucleons would 
contribute to the strange condensate and pure strange 
hadrons would contribute to light quark condensate. 
This effect is captured by the parametric mass depen- 
dence. 
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FIG. 3: Light sigma terms calculated with the constituent 
quark picture. 
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FIG. 4: Strange sigma terms calculated with the constituent 
quark picture. 



To compare with the lattice results of the Wuppertal- 
Budapest group Q the quantity 



A,, a (T) = 



<qq> 



< SS > 



< 11 >o -^r < ss > 



(31) 



is considered. The reason to define this quantity on 
the lattice is purely technical: in this form it elimi- 
nates quadratic singularity at nonzero value of quark 
mass m q /a 2 (where a is lattice spacing) and the ratio 
eliminates multiplicative ambiguities in the definition of 
condensates. Physically this quantity is sensible to chi- 
ral symmetry restoration: it is normalized to unity in 
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FIG. 5: Values of the Ah coefficient for hadrons of different 
mass from formulaf2Z]as evaluated with the generalized quark 
counting formula. 



vacuum and vanishes with vanishing of the condensates 
as temperature grows. Fig. [5] shows a comparison of 
the lattice data to the HRG results with CQP mass 
formulas. The band set by the LO chiral masses for 
pscudo-Goldstone bosons (n = 0) and NLO corrections 
(k = 0.021) encapsulates the 1QCD data. There is an 
overall agreement up to temperatures ~ 155 MeV which 
is the critical temperature from the lattice data. To com- 
pare in Fig. [9] we show HRGM results with parametric 
mass formulas. Already there is less good agreement only 
in temperatures up to ~ 140 MeV. 



VI. HOLOGRAPHIC MASS FORMULAE 

The second model considered in this pap er is the holo- 
graphic model of Sakai and Sugimoto [19|. This model 
is based on a D-brane construction in string theory and 
assumes both large N and large 't Hooft coupling g 2 N. 
Even though this model is neither supersymmetric nor 
conformal, the approximations used are sufficiently un- 
der control to justify the serious effort that has gone into 
exploring its phenomenology. Even though in its orig- 
inal formulation the model did not allow for non-zero 
quark masses, it leads to a large number of quantitative 
predictions which agree very well with experiment de- 
spite the model having just two parameters [48j]. The 
inclusion of explicit chiral symmetry breaking by non- 
vanishing quark masses was studied 3 in subsequent work 



[22], |49j ■ The resulting hadron mass shifts were calculated 
for the case of two flavours in [24j and for three flavours 
in [25|. The latter reference provides hadron mass for- 
mulae which were used in the present study. The results 
reported here include only the nucleon octet and delta 
decuplet states in the sums over hadrons, since mass for- 
mulae have only been calculated for these states. 

In the quasi-Goldstone boson sector the holographic 
model leads to the GMOR formula [22| . Although in 
the holographic model the GMOR relations were only 
obtained in the leading order m 2 = 2c/ f 2 (m u + m,d) and 
vr? K = 2c/ f 2 (m u +m s ) [22j , in the following the equations 
1 181 191 will be used, which include an estimate of higher 
order corrections parameterized in therms of the constant 
K. 

For the baryon sector the results are as follows. In 
the case of two quark flavors the formula for the nucleon 
octet and delta decuplet reads p4| 



cm, 



(32) 



where c = 4.1 GcV~ . The leading order chiral per- 
turbation theory result is of exactly the same form with 
c = 3.6 GeV~ : [44[ (and references therein). For the 
choice of parameters made in this paper this mass shift 
gives 8M = 80.36 MeV and the resulting sigma term is 



a = —cm n 



< qq >o 
f 2 



2k- f ) = 36.8641 MeV . (33) 

J7T 



For an alternative approach see [23 



Note that the above results are state independent. 

The estimated pion-nuclcon sigma term in chiral per- 
turbation theory changes from a^N ~ 59 ± 19 MeV at 
NLO to a wN = 43 ± 7 MeV at N 3 LO which means 
is quite close to what was obtained above (within er- 
ror bars). There is no essential difference for this sigma 
term when one includes the strange quark, which is why 
one can compare this with the 2 + 1 results. In the chi- 
ral limit the nucleon octet baryons was found Q to be 
M = 767 MeV which gives SM = 171 MeV from the 
physical proton mass. 

For the two flavour DSE studies [l2[ the nucleon and 
delta sigma terms were found to be — 60 MeV and 
cta — 50 MeV. This is within the resonable limits defined 
by various model approaches but will turn out to be a 
little closer to the holographic results of 2 + 1 flavour 
case. 

In the three flavour case [25[ the nucleon octet mass 
formula reads 

SM N = ^c 8 (a m 2 K0 + a K m 2 K± + a^m 2 ^) , (34) 

and for the delta decuplet 

SMa = -^c w (a m 2 K o + aKm 2 K± + a v m\±) , (35) 

where c 8 = 7.9 GcV -1 , Cio = 9.5 GeV -1 and the a coeffi- 
cients are given in tables [TTI and IIIII Using equations (f2"0")) 



9 




FIG. 6: HRG model results. Left panel: temperature dependence of the quantity A av defined by the equation (|29|) . Right 
panel: temperature dependence of the thermal reduced mass m rc( j defined by Eq. (|30p . The parameter mo denotes the upper 
limit for the mass of hadrons included in the calculation. 




250 



FIG. 7: HRG model results. Left panel: Quark counting mass formulae based on quark content of a hadron. Green solid 
line light quark condensate green dashed line - strange quark condensate. Right panel: Chiral condensate with parametric 
dependence of hadron masses [27T429I ] . 



and (|22p one can calculate derivatives of baryon masses 3(SMb) 1 



d(SM B ) 



r* 



CLK- 



< qq > + < ss > 



2/2 



K 



1 + 2k 



< qq >o , 

a n — — I 1 



2k- 



p 

■> K 



_K_ 

■I 7T 



(36) 



dm d 



d(SM E 



r* 



< qq > + < ss > / m z K 

ao — m — v +2k it. 



< qq >o (, 



(37) 



dm s 

where B — N, A and # 



1 < qq > + < ss > / m\ 
-c#{a Q +a K ) ^ ^1 + 2k— 



(38) 



10. 
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CQPk=0.021 

• Lattice 
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200 



stituent quark picture. 



The resulting hadronic sigma terms are presented in 
the tables ITVW1 

In the holographic setup the strange nucleon sigma 
term is significantly overestimated, indicating that higher 
order corrections are needed. For ChPT the leading order 
tree level result is expressed in terms of five low-energy 
constants and gives a reasonably good evaluation of the 
nucleon strange sigma term. For the purpose of compar- 
ison, the results for the leading order ChPT sigma terms 
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ao 
an 



P N A E 4 



3/5 4/5 9/10 3/5 11/10 8/5 4/5 8/5 
4/5 3/5 9/10 8/5 11/10 3/5 8/5 4/5 
8/5 8/5 6/5 4/5 4/5 4/5 3/5 3/5 



TABLE II: Coefficients in the nucleon mass formula 
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A++ 


A+ 


A 


A" 


£*+ 


E*° 


£*- 






Q~ 


to the lattice 




1/2 


3/4 


1 


5/4 


3/4 


1 


5/4 


1 


5/4 


5/4 


CQP - con- 


a K 


5/4 


1 


3/4 


1/2 


5/4 


1 


3/4 


5/4 


1 


5/4 




a n 


5/4 


5/4 


5/4 


5/4 


1 


1 


1 


3/4 


3/4 


1/2 




TABLE III: Coefficients in the delta mass term 



150 200 
T [MeV] 

FIG. 9: Comparison of the HRG results to the lattice 
Wuppertal-Budapest results [3] (blue dots). Red, green lines 
parametric mass formulae. 



are shown in table HI The strange sigma term for the 
nucleon, er^ w 162 MeV, is a bit large but still reason- 
able. It should be noted that when compared to the NLO 
results from [14j even the sign of the sigma terms can 
change, meaning that higher order corrections cannot be 
ignored. It is to be expected that including higher order 
corrections in the holographic approach should cure the 
problem of overestimating the strangeness contribution 
as it does in the case of ChPT. 

Fig. [10] presents the result for the chiral condensates 
obtained with holographic and NLO ChPT mass for- 
mulae where only the nucleon octet and delta dccuplet 
baryons are included (apart from the quasi- Goldstone 
bosons). In the holographic case, due to the overesti- 
mated sigma terms, the difference between strange and 
light condensates is diminished. For the same reason the 
too small number of states included in the strange sector 
is compensated. 
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P 


N 


A 


£+ 


£° 


E - 






a s 


50.5611 
47.3916 
1189.1 


47.3916 
50.5611 
1268.62 


42.6751 
42.6751 
1070.76 


44.2976 
28.4501 
713.838 


36.3738 
36.3738 
912.652 


28.4501 
44.2976 
1111.47 


39.5622 
26.8842 
674.548 


26.8842 
39.5622 
992.651 



TABLE IV: Sigma terms for the nucleon octet in MeV. 
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A+ 


A A- 
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E*~ 






fT 


&d 
<J a 


49.4056 
37.5199 
941.409 


45.4437 
41.4818 
1040.82 


41.4818 37.5199 
45.4437 49.4056 
1140.22 1239.63 


43.4863 
35.5626 
892.297 


39.5245 
39.5245 
991.705 


35.5626 
43.4863 
1091.11 


37.5671 
33.6052 
843.186 


33.6052 
37.5671 
942.593 


31.6479 
31.6479 
794.074 



TABLE V: Sigma terms for the delta decuplet in MeV. 




FIG. 10: Temperature dependence of strange and light condensates using holographic mass formulae (left panel) and NLO 
ChPT mass formulae [l4| (right panel). For details, see the text. 



As is clear from the above discussion, one important 
extension of the existing calculations in the holographic 
model would be to calculate higher order corrections in 
the current quark masses. One motivation for it was al- 
ready mentioned: this would improve the resulting sigma 
terms, especially in the strange sector. A second, more 
formal, motivation is that in ChPT for Nt — 2 the sec- 
ond order correction to the proton mass has the universal 
form [|U 

= M + 4 Cl ml + , (39) 

where §a is the axial coupling. If one adopts the usual 
scaling of parameters with TV, then one gets qa ~ N and 
f%^N,so that the sublcading contribution would scale 
like ~ N 2 which would dominate the leading order re- 
sult Mq ~ TV. On the other hand, if one follows recent 
argumentation [5(j that <?a ~ TV = 1 then NLO con- 
tributions would be of the order ~ \/N. It would be 
interesting to check if in the Sakai-Sugimoto model this 
universality also holds. 



VII. CONCLUSIONS 

This paper was devoted to a discussion of the finite 
temperature behaviour of the chiral condensate within 
the HRG framework exploring different microscopic de- 
scriptions of the dependence of hadron masses on the 
current quark mass. In particular, a constituent quark 
scheme and holographic mass formulae have been used. 
It was also studied how the results are affected by in- 
cluding different numbers of states in the sums over res- 
onances. It turns out that with a sensible choice of mass 
formulae and including hadron states with masses up to 
~ 2 GeV generic agreement with recent lattice results is 
obtained. This is yet another confirmation of the well 
known fact that for low temperatures the HRG model 
gives a satisfactory physical interpretation of 1QCD data. 
Chiral symmetry restoration in the strange sector was 
seen to take place at higher temperatures than in the light 
quark sector, which is related both to the lower number 
of strange hadrons contributing to the condensate as well 
as to the response of hadron masses to changes in the cur- 



rent strange quark mass. 8975/E-344/M/2012 for young scientists. 

A generalization of the quark-counting approach of 
[Til fig was proposed, and it was shown that the mass 
relations where only valence quarks of the hadron are 
taken into account already lead to a behaviour of the con- 
densate which is close to what is seen in the full lattice 
data. In this scheme dynamically generated (constituent) 
quark masses are considered and their dependence on the 
current quark mass is quantified in the framework of the 
NJL model. This step takes into account part of the 
non-perturbativc QCD dynamics. The sea quark con- 
tributions are neglected resulting in a vanishing strange 
sigma term for the nucleon and a vanishing light quark 
contribution for the f2~ baryon. This is somewhat in the 
spirit of the large- N expansion where quark loops are 
suppressed. 

Along with this, a careful analysis of the hidden 
strange mesons has been performed based on the flavour 
symmetry structure of the mesons. This affects the sim- 
ple quark counting rules used by [HI, HU, taking into 
account neglected effects which overestimated the light 
quark condensate and underestimated the strange quark 
condensate. 

Another new aspect considered in this paper concerns 
the sigma terms and the condensate following from the 
mass formulae of the holographic model of QCD due to 
Sakai and Sugimoto [l9| . These formulae take on a form 
similar to the tree level ChPT results with strange sigma 
terms overestimated due to the inaccuracy of the approx- 
imation for the relatively large value of m s . Since those 
shifts were only calculated for the nucleon octet and delta 
decuplet baryons, the resulting condensate is incomplete. 
This also shows the importance of heavier hadrons for 
temperatures near the QCD transition temperature. 

The results obtained here are of great importance in 
the context of hadron production under extreme condi- 
tions in heavy-ion collisions. Recently, it has been con- 
jectured that the behaviour of the chiral condensate de- 
termines the collision rates of hadrons and thus may pro- 
vide a microscopic approach to the chemical freeze-out 
of hadron species [16j . This approach, however, has yet 
been considered only in the light quark sector. Including 
the strange quark condensate to that model could ad- 
vance the understanding of strangess production in heavy 
ion collision experiments. 
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